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The stability problem of the rotation induced electrostatic wave in extragalactic jets is presented. 
Solving a set of equations describing dynamics of a relativistic plasma flow of AGN jets, an expression 
of the instability rate has been derived and analyzed for typical values of AGNs. The growth rate was 
studied versus the wave length and the inclination angle and it has been found that the instability 
process is much efficient with respect to the accretion disk evolution, indicating high efficiency of 
the instability. 
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I. INTRODUCTION 

Considering the problem of the escaping radiation from 
active galactic nuclei (AGN), one has to note that we 
have a strong observational evidence of a complex picture 
of emission spectra, which starts from the radio, the opti- 
cal up to X-ray and 7-ray, with the bolometric luminosity 
power of AGNs lying in the range: ^ [W^^ — lO^^Jerg/s 
(see Ref. 1). Origin of emission is supposed to be energy 
of an accreted matter and a spinning black hole, however, 
the details of the conversion process of this energy into 
electromagnetic radiation is still unknown. The discov- 
ery of blazars (see Ref. 2) as sources of ultra-high energy 
radiation, has revealed that the formation and accelera- 
tion of relativistic jets are key processes in understanding 
the energy conversion mechanism. 

An innermost region of AGNs is characterized by rota- 
tional motion, and it is obvious that the rotation affects 
acceleration of plasmas in jets, which consequently may 
influence a process of radiation. The problem of a role 
of centrifugal acceleration (CA) as the rotational effect, 
on a theoretical level has been studied in Ref. 3 for the 
Schwarzschild black hole and it has been found that un- 
der certain conditions the centrifugal force may change its 
direction from centrifugal to centripetal. The same effect 
was shown in Ref. 4 where authors studied dynamics of 
particles, sliding along rotating straight channels. It has 
been argued that due to the centrifugal force, the process 
of acceleration might be very efficient for rotation energy 
extraction. On the other hand it is obvious that no phys- 
ical system can preserve rigid rotation nearby the light 
cylinder surface (LCS) (an area where the linear velocity 
of rotation is equal to the speed of light) where mag- 
netic field lines must deviate from the straight configu- 
ration. Generalization of Ref. 4 for curved field lines has 
been presented in Ref. 5 and it was found that for the 
Archimedes spiral a centrifugal outflow may reach infin- 
ity, avoiding the light cylinder problem. According to the 
spin paradigm (stating that, to first order, it is the nor- 
malized black hole angular momentum, that determines 
whether or not a strong radio jet is produced) (see Ref. 
6) rotation is supposed to be fundamental in formation 
of outflows. Blandford & Payne (see Ref. 7) consider- 



ing the problem of creation of the outflows, pointed out 
that in this process the centrifugal acceleration may play 
an important role depending on the inclination angle of 
magnetic field lines. A series of works dedicated to effects 
of rotation, examine its role in the nonthermal emission. 
One of the important examples from this point of view is 
the work of Blandford & Znajek (see Ref. 8), where they 
show that the rotational energy of the black hole can be 
extracted electromagnetically as a Poynting flux. In the 
context of AGN, Gangadhara & Lesch (see Ref. 9) con- 
sidered the problem for understanding how efficient the 
rotation is for producing X - ray and 7 - ray photons. 
It was shown that nearby the light cylinder surface, cen- 
trifugally accelerated particles, due to scattering against 
low energetic photons, may provide the mentioned ra- 
diation. Reconsidering the same problem for explaining 
the TeV energy emission from a certain class of blazars, 
in Ref. 10 it has been argued that the CA might be so 
strong, that could provide ultra-relativistic electrons. 

All mentioned cases exhibit the possibility of energy 
pumping from rotation and in cases of the centrifugal ac- 
celeration they show its high efficiency in the energetics 
of relativistic electrons in rotating plasmas accelerating 
up to energies corresponding to the Lorentz factors of the 
order of 10^~^ (see Ref. 10). Thus, the amount of energy 
contained within a centrifugally accelerated plasma flow 
is very big and if there are mechanisms for the conver- 
sion of at least a small fraction of this energy into the the 
variety of instabilities - one might find a number of in- 
teresting consequences related to the problem of plasma 
energy conversion into radiation. 

Generally speaking, for the energy pumping process, 
one has three fundamental stages: (a) energy conversion 
from rotation into kinetic energy of the plasma back- 
ground flow, (b) kinetic energy conversion into energy 
of electrostatic waves and (c) conversion of the energy 
of waves into radiation by means of the process of the 
inverse Compton scattering (ICS) of electrostatic waves 
to photons I + e = I' + t, which is thought to be re- 
sponsible for creation of nonthermal radiation (see Ref. 
11). By Tautz & Lerche (see Ref. 12) the electro- 
static/electromagnetic mixed mode has been considered 
for relativistic jets in the context of Gamma-Ray Bursts, 
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to determine low-frequency instabilities. But our aim is 
to consider only electrostatic waves and investigate the 
second stage of the process [see (b)] and show how effi- 
cient it is. Therefore the electrostatic wave, playing an 
important role in the process of radiation might be very 
effective if one makes it unstable and as we will see the 
centrifugal force may be very sufficient in this context. 

The centrifugally driven parametric instability has 
been introduced in Ref. 13 where it was shown that 
the centrifugal force acting on particles inside the pulsar 
magnetospheres, causes separation of charges, leading to 
generation of an electrostatic field, which excites the cor- 
responding instability. The increment of the linear stage 
has been estimated and analyzed for the Crab pulsar and 
it turned out that the linear regime was extremely effi- 
cient and short in time, indicating a need of saturation of 
a growth rate. This mechanism is called parametric, be- 
cause an "external" force - the centrifugal force, playing a 
role of the parameter, acts on plasma particles, changes 
in time and gives rise to the instability. By applying 
the same approach, in Ref. 14 (to be published in Mon. 
Not. R. Astron. Soc.) the problem of the reconstruction 
of the pulsar magnetospheres on large scales nearby the 
LCS has been considered. We have studied a new para- 
metric mechanism and it was shown that for curvature 
drift waves, corresponding instability might be extremely 
efficient. 

In this paper a hydrodynamic approach will be used 
to study the parametric mechanism of rotation induced 
electrostatic wave generation in extragalactic jets, apply- 
ing the approach developed in Ref. 13 and Ref. 14 (to 
be published in Mon. Not. R. Astron. Soc). 

The paper is arranged as follows. In SjlT]we derive the 
dispersion relation, in mill the corresponding results are 
present and in ^ilVI we summarize our results. 



II. THEORY 



tion of motion of a particle is given by: 

^=7g+-(E + vxB). (2) 
ar m 

where 7 = (1— v^)~^/^ is the Lorentz-factor, v — dr/dr is 
the velocity of the particle determined in the 1-1-1 formal- 
ism, p — > p/m - the dimensionless momentum, g = — 
7g - the effective centrifugal force (see Ref. 13) and 
T = (^ = ^1 — ilg^-r^) - the universal time. Tak- 
ing into account an identity d/dr = d/{£,dt) -I- (v- V) and 
a fact that Lorentz factors in the co-rotating frame of ref- 
erence and in the inertial-laboratory frame (LF) relate to 
each other: 7 — ^7' (prime denotes a physical quantity 
in the LF) , one can rewrite the equation of motion in the 
LF by following: 

^ + (vi-V)pi = -7eV^+-(E + VixB), (3) 
ot m 

complemented by the continuity equation: 

^ + V-(n,vO=0, (4) 
and the Poisson equation: 

V • E = 47r^e,nj. (5) 



i — e,p,b 

here e,p, b denote electrons, positrons and the bulk com- 
ponents respectively. Eq. ([3]) in the zeroth approxima- 
tion (leading state), by applying the frozen- in condition 
Eo + voi X Bq = 0, can be reduced to (see Ref. 4): 



Throughout the paper it is supposed that magnetic 
field lines are straight and inclined by the angle a with 
respect to the rotation axis. The plasma is in the frozen- 
in condition, co-rotating together with the field lines with 
the angular velocity fi. 

It is easy to start our consideration by the local non in- 
ertial co-rotating frame of reference. The corresponding 
interval in the rigidly rotating frame will have the form 
(see Ref. 10): 



1 - Vl 



(6) 



In Ref. 4 it has been shown that for ultra-relativistic 
particles (7^1) Eq. ((6]) has following solutions: 



r{t) = sin^e ft, 



(7) 



ds^ = - {1 - nlf-r^) - dr^ , (1) 

where ile/ = Q sin a is the effective angular velocity of 
rotation. Here we use units c = 1. In this paper we apply 
a method developed in Ref. 13 where relativistic plasma 
motion was described by the single particle kinematics. 
According to this approach, by using 1-f 1 formalism (see 
Ref. 15), in the co-rotating frame of reference the equa- 



vo{t) = Vbcosr^e/i, (8) 

for initial conditions: rg = and Vq ~ 1. As we see from 
Eq. (I3|), the effective centrifugal force — 7^^^ changes 
in time due to Eqs. (|7l8p and therefore the parametric 
mechanism switches on, giving rise to the electrostatic 
instability (as we will see). 

When studying the problem for the leading terms, 
plasma oscillations have not been considered and only 
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the effects of the centrifugal acceleration have been taken 
into account (see Ref. 4). Generally speaking, differ- 
ent species of plasmas [electrons, positrons and protons 
(bulk)] experience the centrifugal force, as a result they 
separate, which will lead to the generation of the electro- 
static field, creating the Langmuir waves. This process 
can be considered as a next step of the approximation. 
The aim is to investigate the linear perturbation theory 
of the electrostatic instability and estimate its role in rel- 
ativistic AGN jets. 

We start the analysis by introducing small perturba- 
tions around the leading state: 



70£ 



-^y^Js{a)Nb{io + s^ef), 



(18) 

where M and m are masses of protons and elec- 
trons/positrons respectively, a = ck/VLef and = 
y'SvmlJe^/TO - electron/positron plasma frequency. For 
deriving Eqs. (|17ll8p . the identity 



^zEix sin fieji 



(19) 



where ^ = (n,v,p,E,B). 

Perturbing all physical quantities by following: 



*^(t,r) cx *^(i)exp[i (kr)] , 
Eqs. will get the form: 



(9) 



+ ikvon] + ikuiovl = 0, 



ikE^ = in'^etnl, 



(10) 

(11) 
(12) 
(13) 



Here is the electric field in the LF, induced by sepa- 
ration of charges. 

In order to reduce the above system into a single equa- 
tion let us use an ansatz: 



n} = N,e 



then from Eqs. (|11I12I13P one gets 



at2 



(14) 



(15) 



has been used. Direct substitution of Eq. (fT7|) into Eq. 
(Ull), yields: 



^6 



Jsia)Ji{a) 



LU + (s — l)^ef 



X N{lO + (s - l)Qef)- 



(20) 



where Wf, = ^/87^rl^e^JM is the plasma frequency corre- 
sponding to the bulk component. 

Generally speaking in order to solve Eq. (PO)) one 
has to consider similar equations, rewriting Eq. ([^0]) for 
N{u!±ilef), A^(w±2f2e/), etc., thus for solving the prob- 
lem exactly, one needs to solve a system with the infinite 
number of equations, which makes the problem impossi- 
ble to handle. Therefore the only way to overcome this 
difficulty and to gain extended view of the general be- 
havior of the instability, is to consider physics close to 
the resonance condition, which provides the cutoff of the 
infinite row in Eg. ((20)) and makes the problem solvable 
(see Ref. 16). 

If one considers in Eq. (|20p only resonance terms, cor- 

1 /2 

responding to the following frequency tJres ~ <^e/7oe j 
then preserving leading terms and taking into account the 

resonance conditions: uj w —so^ef, sq = Iq 



(here [A] means the integer part of A), the dispersion 
relation can be reduced into a following single term spec- 
ifying the growth rate (A ) of the instability: 



ikE^ = 47r^e,iV,e-*-^% 



(16) 



where = ^ sin(f]e/t). 

Restoring the speed of light and introducing a new 
variable N = Np — N^, after making the Fourier trans- 
formations, one can easily reduce Eqs. ()15|16|) : 



(17) 



1/2 '^so\"')- 



47ofc7o, 



(21) 



Since we are interested in imaginary parts of A, it is 
easy to see that the following solution (see Ref. 17): 



1/2'^ so\"') 



47o67o 



(22) 



is responsible for the instability, the increment of which 
is given by: 
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FIG. 1: Dependence of logarithm of the instabihty rate on 
the wave length. The set of parameters is: a = 1°, Mbh ~ 
W^Mq, n = 3x W-^s-\ Job = 20, 7oe = 10^ and = n° = 
O.OOlcm"^. As it is seen in the figure, the increment is very 
sensitive on the wave length: slightly changing A, one may 
kill the instability completely. 
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III. DISCUSSION 
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FIG. 2: Dependence of logarithm of the average instability 
rate on the wave length. The set of parameters is the same as 
in Fig. [T] except a wider range of a and fixed A = Ao = lO^cm. 



their ambient, obviously having density of the order of 
n-am ~ lcm"3. 

An interesting feature of the result shown in Fig. [1] 
is sensitiveness of the growth rate on the wave length. 
As it is clear from the plot, a small change of A, drasti- 
cally changes the situation: for certain values of the wave 
length the increment reaches its maximum level and for 
slightly different values - it becomes equal to 0. Therefore 
it is better to examine an average value for each interval 
with a single peak, taking into account both parameters: 
A and a, calculate the growth rate with respect to them: 



We consider a central black hole with mass Mbh = 
IO^Mq and an angular rate of rotation = 3 x lO^^s^^ 
of an AGN wind, making the light cylinder located at a 
distance ~ lO^^cm from the center of rotation. These 
values are typical for active galactic nuclei (see Ref. 18). 
AGN winds and extragalactic jets are supposed to be 
composed of a bulk component having the Lorentz factor 
of the order of ^ 20 and relativistic electron-positron 
pairs with Lorentz factors: ~ 10^ (see Ref. 10). 

As it was already explained, the centrifugal force leads 
to the generation of the Langmuir waves with the wave 
vector k {\\B), the growth rate of which we are going 
to estimate versus the wave length A and the inclination 
angle a in order to investigate the corresponding insta- 
bility. 

Let us suppose a jet with the opening angle (3 = 2a = 
2°. In Fig. [1] we show the dependence of logarithm of 
the instability increment on the wave length A = Ao + 
AA, where A = lO^cm and AA = [0,0.1]cto. The set 
of parameters is: a = 1°, Mbh = 10® Mq, f7 = 3 x 
10-^-1, 7ob = 20, 7oe = 10^ and n° = n° = O.OOlcm^^. 
The values of n° and rig are typical for extragalactic jets, 
which are supposed to be under dense with respect to 



r ^ n VT( V r r d^daSiX, a), (24) 

(A2 - Ai)(a2 - ai) J^^ J^^ 

and based on discrete data interpolate it for a wider range 
of parameters. Here A1.2 and ai^2 are minimum and max- 
imum values of the wave length and the angle respectively 
for each interval. 

In spite of that jets usually have small opening angles, 
it is better to study a dependence of the growth rate on 
inclination for a wider range of angles, in order to under- 
stand a general behaviour of the increment. In Fig. [2] 
we show the dependence of logarithm of the growth rate 
versus a. The set of parameters is the same as for Fig. [T] 
except a wider range of a and a fixed value of the wave 
length A = lO^cm. As we see, the increment is a con- 
tinuously increasing function. This is a natural result, 
because when one increases the angle, the effective angu- 
lar velocity increases as well, that makes the centrifugal 
force higher, leading to the more efficient instability pro- 
cess. 

As we see from the present figure, the instability is 
less efficient for smaller angles. Let us study now a de- 
pendence of the growth rate on the wave length for the 
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FIG. 3: Dependence of logarithm of the average instabihty 
rate on the wave length. The set of parameters is the same 
as in Fig. [TJ except a wider range of the wave length. 
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FIG. 4: Dependence of logarithm of the evolution time scale 
on the luminosity power. The set of parameters is: r; = 0.03, 
e = 0.1 and Mbh = 10* Mq. luminosity power lies in the 
range L e 1.4 x [10**; 10'"^]erc//s. 



inclination angle 1° when the instability is minimal (as 
it it is seen in Fig. [2]). 

In Fig. [3] the behaviour of logarithm of F versus A is 
shown. The set of parameters is the same as in Fig. [T] 
except a wider range of the wave length. The figure shows 
that by increasing A, the average growth rate increases 
as well. Let us consider a less efficient case, thus the 
smallest increment (for A = lO^cm) shown in the figure: 
F ^ 1.3xl0~^s~^ (A = lO^cm), and see how powerful the 
instability is. One can easily estimate the corresponding 
time scale by inversing the value of F: 



(25) 



On the other hand jets are formed due to an accret- 
ing matter, therefore energetics of jets is defined by the 
efficiency of the accretion process, and in order to under- 
stand how efficient the considered instability is, one has 
to estimate the evolution time scale of an accreting disk 
and compare it with the instability time scale. 

In Ref. 19 the problem of fueling of AGNs is consid- 
ered and it is shown that for a given mass and luminosity 
of the AGN, mass of the self gravitating disk can be ex- 
pressed by: 



to the accretion mass rate M and the luminosity L by 
the following expression: 



(27) 



Mc2 

Lb is the Eddington luminosity 



Le = 1.4 X 10« J|^er5/s. 



(28) 



Then estimating the disk's evolution time scale tevoi = 
Mgg/M, by taking Eq. into account, one can easily 
get: 



= 3.53x10 
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(29) 



Msg = 2.76 X 10^ 
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(26) 



where ry is the standard Shakura, Sunyaev viscosity pa- 
rameter (see Ref. 20), e is the accretion parameter related 



For taking the sense out of the considered parametric 
instability, let us examine typical values of the accretion 
disk 7] = 0.03, e = 0.1 and Mbh = 10^ Mq, then consid- 
ering the dependence of the evolution time scale on the 
luminosity, one gets the plot shown in Fig. |4l For produc- 
ing the plot we have examined the following luminosity 
range L/Le & [0.01; 1]. As we see in the figure, the evo- 
lution time scale is a continuously decreasing function, 
and for the minimum value of the considered luminosity, 
the time scale is of the order of lO^^s, whereas for the 
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Eddington limit the time scale reaches its minimum level 
~ 10"s. 

Calculating tevoi/tinst, by using Eq. ([25| one can see 
that the present ratio lies in the range ~ [10^; lO^^Js, 
which means that the process of the instability is much 
faster than typical rates of accretion, therefore the pro- 
cess of conversion of rotational energy into energy of the 
Langmuir waves is very rapid. This means that the linear 
stage has to be short in time, and the growth rate must 
saturate soon due to non linear effects, consideration of 
which seems to be essential. 

In the introduction we have noticed that the Lang- 
muir waves may strongly influence processes of radia- 
tion through the ICS of these waves against soft photons, 
when due to the channel I + e = I' + t the electrostatic 
waves lose energy, whereas the photons gain energy. This 
problem was not an objective of the present paper, which 
seems to be the first step for understanding the processes 
of rotation energy pumping into instabilities. As a next 
step one has to study this particular problem as well and 
see how efficient the radiation (produced via the ICS) 
could be for the extremely unstable electrostatic waves. 
This will be an objective of a work which we are going 
to make soon or later. 



IV. SUMMARY 

1. We studied a relativistic plasma of extragalactic 
jets composed of the bulk flow and relativistic 
electron-positron components. The role of the cen- 
trifugal acceleration in generation of the parametric 
instability of the Langmuir waves has been consid- 
ered. 

2. Perturbing the Euler, continuity and induction 
equations, preserving only first order terms, we 
have derived the dispersion relation governing the 
instability. 

3. Examining the resonance condition of the system, 
an expression of the growth rate has been obtained 
for a nearby zone of the light cylinder surface. 

4. Studying dependence of increment on the wave 
length for typical values of extragalactic jets, we 
found that the instability becomes efficient only for 
certain and narrow ranges of the wave length. 



5. Taking into account efficiency of an accretion pro- 
cess in AGNs, and comparing the corresponding 
time scale to the time scale of the electrostatic in- 
stability it has been found that even when the latter 
is less effective (small angles and wave lengths) the 
instability process is extremely efficient. 

The electrostatic instability may strongly affect 
processes of radiation in AGN jets by means of the 
ICS, therefore for extending this work, it is natural 
to consider a next important step and study the 
last stage: energy conversion into radiation. 

An important restriction in the present model is 
that we studied the problem for quasi-straight mag- 
netic field lines, whereas in real astrophysical situ- 
ations the field lines are curved and it is interesting 
to see what changes in the dynamics of the electro- 
static wave instability, when the curvature is taken 
into account. 

On the other hand the curvature of magnetic field 
lines induces the curvature drift waves, which also 
may affect the process of energy pumping. In the 
context of pulsars this problem has been consid- 
ered in Ref. 14 (to be published in Mon. Not. 
R. Astron. Soc). It is reasonable to examine the 
same problem for relativistic AGN jets and study 
the role of the curvature drift modes in the energy 
conversion process. 

When studying the problem on the linear stage, the 
increment shows high efficiency of the electrostatic 
instability, therefore the non linear regime also has 
to be examined in order to study the problem for 
longer time scales when the linear approximation 
does not work any more. For this purpose we plan 
to implement a numerical magnetohydrodynamical 
code to study this particular case as well. 
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